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For a class of quasi-one-dimensional clusters, by using exact diagonalization, we study the effect 
of side spins on the spin-spin correlations on chain. Our calculations show that the side spins added 
in the same sublattice can effectively strengthen the spin-spin correlations in large distance region 
and make the change tend to flat. It is exactly proved that periodically adding side spins can set 

■ up magnetic long-range orders in the ground state. Also we investigate the effect of the density of 
side spins on correlation strength. The case that two sublattices have different localized spins is 
discussed. 

Spin-spin correlation and magnetic long-range order (LRO) are of fundamental importance for quantum spin sys- 
tems. They have been studied by many analytical and numerical works. For the bipartite Heisenberg antiferromagnetic 
(AF) systems, the spin-spin correlations of any two sites in the ground state (GS) always are AF, i.e. the correla- 
tion functions are positive when two sites belong to the same sublattice whereas negative the different sublattices. 
However, it is not sufficient to set up the AF LRO. The dimensionality of system is one of the key factors. It has 
been rigorously proved that there exists Neel order in the GS [Q of the three-dimensional system, whereas no Neel 
order for the one-dimensional one. Although there is no magnetic LRO for one- and two-dimensional (ID and 2D) 
Ch ' Heisenberg antiferromagnets at temperature T > by Mermin- Wager theorem T = may be the critical point. 
I | When the localized spin S > 1, it has been proved that the GS of the 2D Heisenberg antiferromagnet has Neel order 
H . Notwithstanding the rigorous proof of Neel order when S = \ has not been established until now, many numerical 
and analytical works support its existence at T = Q| . 

The situation of quasi-one-dimensional (QOD) systems may be diverse due to their various geometric structures. 
The discovery of a heavy- fermion phenomenon in the Ce-doped Neodymium cuperate [^) has led to an increasing 
interest in the study of strongly correlated electrons coupled antiferromagnetically to magnetic moments . W. 
Zhang et al. have investigated the case of a single magnetic impurity and found that the spin-spin correlation function 
between the impurity spin and spins in the chain extends over long range 0. Obviously, the magnetic LRO can not 

■ be established through adding a single impurity spin to an ID chain. The situation may change when impurity spins 
are periodically added to the chain in a special way. In this paper, by using exact diagonalization, we investigate a 

f— \ , class of QOD Heisenberg AF clusters, in which there are some impurity spins sit beside an ID finite chain (we call 
y—i ' them as side spins), and explore the effects of side spins on spin-spin correlations. Our numerical results indicate 
that the spin-spin correlations in the region of large distances can be enhanced by adding side spins in the same 
sublattice. When periodically adding side spins, the decay of spin-spin correlations with distances slows down and 
becomes obviously flat in the range of large distances. For infinite ID chains, it is analytically proved that adding 
side spins can set up magnetic LROs. Also, we investigate the variation of magnetic LROs with the density of side 
spins, and find that the decay of ferromagnetic (F) LRO is faster than that of AF LRO. 
. We investigate the spin-i Heisenberg AF system with interactions of nearest neighbors, whose Hamiltonian reads 
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H = jJ2S(i)-S(j), (1) 
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where J > 0. (i,j) denotes the sum over pairs of nearest neighbors. The spin-spin correlation can be written as 

A(Ri — Rj) = (g S(i) ■ S(j) GJ. Here, Ri and Rj are the coordinates of sites i and j respectively, and \G) represents 

the GS. Let i denote the site by the ith site of the chain. At first, we consider a finite chain of 23 sites under free 
boundary conditions and study the effect of a single side spin on the spin-spin correlations. By exact diagonalization, 
we calculate A(R t — Rj) for site 2 (i.e. j = 2 and i = 3,4, 5, ...), and the numerical results are shown in Fig. 2, here 
the side spin is placed at site 2, 3 and 4 respectively (see Fig l.(a)). When the side spin is added on site 3, A(Ri — Rj) 
becomes weaker than that of the pure ID chain. Namely, the side spin weakens the spin-spin correlations. But, when 
the side spin is at site 2 (or 4), A(Ri — Rj) becomes stronger in the region of large distances. In this case, the side 
spin strengthens the spin-spin correlations for most of (Ri — Rj). When a single side spin is placed at sites 5, 7, 9 
and 11 respectively, our calculations show that the effect of the side spin is similar to that of the side spin at site 3. 
And, when the side spin is placed at site 6, 8, 10 and 12 respectively, it is similar to that of the side spin at site 2. 
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Under the periodic boundary conditions, the case of a single side spin has been studied by Monte Carlo simulations 
and the spin correlation function is found to extend over long range . For such a system, the GS always takes the 
lowest possible total spin (LTS) no matter where the side spin is attached. But the situation for a finite chain under 
free boundary conditions becomes little complicated. The finite chain can be divided into two 'sublattices' as doing 
for infinite systems. Sites 1, 3, 5, 21 and 23 belong to sublattice A, and sites 2, 4, 6, 20 and 22 sublattice B 
(Fig. 1(a)). Adding the single side spin at different sites can leads to the GS taking different values. When the side 
spin sits by the site of sublattices B 7 it belongs to sublattice A. By Lieb-Mattis theorem pQj , the GS has the global 
spin 1 and is 3-fold degenerate. In other words, the GS takes a higher total spin than its lowest possible value. But, 
when the side spin sits by the site on sublattices A, it belongs to sublattice B. The GS takes the LTS zero. Our 
calculations show that the effect of a single side spin depends on where the side spin is. In other words, the total 
spin of the GS is related to the behaviors of spin-spin correlations. One important question concerned is how the 
spin correlations change when more side spins are added. When the side spins added are in the same sublattice, the 
total spin of the GS becomes higher, and the difference between it and the lowest possible total spin larger. From the 
above calculations, we speculate that A(Ri — Rj) will become stronger in the region of large distances. 

For further details, we consider the cases of two side spins. If adding two side spins at site 2 and 4, the GS has the 
global spin | while the lowest possible total spin is |. But, if adding two side spins by site 2 and 3, the GS takes the 
lowest possible total spin ~ . From the above notion, one can speculate that spin-spin correlations will become stronger 
if adding two side spin at sites 2 and 4, whereas those weaker at sites 2 and 3. We calculate the spin-spin correlations 
between site 2 and others. The numerical results are shown in Fig. 2. They agree with the above speculation. 

Also, we calculate spin-spin correlations for the following three cases (see Fig. 1(a)): 1). 4 side spins sitting at 
sites 2, 8, 14 and 20 on a finite chain of 21 sites; 2). 5 side spins sitting at sites 2, 6, 10, 14 and 18 on a finite 
chain of 19 sites; 3). 8 side spins sitting at sites 2, 4, 6, 8, 10, 12, 14 and 16 on a finite chain of 17 sites. Their 
ground states have the global spins S T — I, 3 and § respectively. The numerical results are plotted in Fig. 3. 
Our calculations show that periodically adding side spins in the same sublattice can obviously enhance the spin- 
spin correlations in the large distance region. And, as the density of side spins increasing, spin-spin correlations 
become stronger and decay more slowly. Approximately, we can fit the spin-spin correlations in exponential way, i.e. 
A(Ri — Rj) = Ccxp(— \Ri — Rj\ /(), here C is coefficient and £ the correlation length. We calculate C by A(i? 9 — R 2 ) 
and A(R 15 -i? 2 ). For the bare chain of 21 sites, C = 9.30. But, for the above three cases, ( = 51.38, 77.27 and 7345.38. 
It shows that the large distance behavior of spin-spin correlations is enhanced by adding side spins. Especially, £ ~ 10 3 
for the third case. It is much larger than the correlation length of bare chain, even the size of system. Consequently, 
side spins can slow down effectively the decay and make the variation become flat in the region of large distances 
(Fig. 3). It seems to exist the AF LRO. 

Now, we turn to consider an infinite QOD system, which is constructed by periodically adding side spins in the 
same sublattice (Fig. 1(b)). The total number of sites is N — K(l + 2), here K denotes the number of cells and I the 
number of sites on chain between every two side spins. I must take odd, i.e. I — 2k + 1, here k = 0, 1, 2.... Supposing 
the side spin is in sublattice A, the number of sites of sublattice A is Na — K(l + 3)/2 and that of sublattice B 
N B = K(l + l)/2. By the Lieb-Mattis theorem |l(J, the global spin of the GS is 

A = K\(k + 2)S A -(k + l)S B \, (2) 

where Sa and Sb are the values of localized spins on sublattice A and B respectively. To investigate the existence of 
magnetic LRO, one needs to calculate the quantity 

0(q) = (G|£(-q)-£(q)|G). (3) 

S(q) = S(j) exp(iq • j), q is a reciprocal vector. The criterion of magnetic LROs is that <?(q) > O(N) at some 

q- If g(Q) > O(N), here Q = (tt, tt, tt), there exists AF LRO. If g(0) > O(N), there exists F LRO. Following the 
approaches developed by G. S. Tian |12|l|, one can obtain 



A 2 + A 

g(Q) >g(0) =— (4) 

We discuss the case that the sites on two sublattices have the equal localized spin, i.e. Sa = Sa = S. From equation 
(2), we readily obtain A = NS/ (1 + 2). As long as I is finite, it is always true that 

g(Q) >g(0) >0(N). (5) 

From these inequalities, one can conclude that AF and F LROs coexist in the GS, and the former is predominant. In 
other words, side spins set up magnetic LROs. 



2 



Obviously, although there always exist magnetic LROs for finite I, both of F and AF correlation strengths will 
depend on the density of side spins, which is defined as r\ = 1 /(I + 2). We introduce the following two quantities to 
measure F and AF correlation strengths respectively, 



G 



(6) 



and 



AF 
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where Ay = 1 when sites i and j belong to the same sublattice and Ay = — 1 when sites i and j the different sublattices. 
Tp and Taf will decrease as / increasing. When I reaches to the infinite, the system changes into an ID chain and the 
magnetic LROs vanish (i.e. Tf = Taf = 0). Possibly, the variation speeds of Tf and Taf are different. We calculate 
the ratio p = Tf/Taf for I = 1, 3 and 5, and plot the data in Fig. 4. p decreases as rj decreasing approximately 
in the linear way. Then the decay speed of F correlation strength is faster than that of AF correlation strength. In 
other words, the ferrimagnetism becomes weaker and weaker as the density of side spins decreasing. 

Ref. (l3) has proved that the GS of ID Heisenberg AF chain has magnetic LROs when its two sublattices have 
unequal localized spins. It is interesting to investigate the QOD system with unequal localized spins. From equation 
(2) and (4), we can conclude that if 
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S b 
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g(0) = due to A = 0. It means that there is no F LRO. The simplest case is Sa = \ , Sb = 1 and I = 1. One 
can give a spin picture of the GS in valence-bond version. The spin on sublattice B can be divided into two ^ spins. 
One of them forms a singlet with the nearest side spin, and the other combines with its nearest neighbor on chain. 
We think this kind of configurations governs the physics of the GS. And it is responsible for the F LRO to vanish. 
Although we have not exactly proved that the AF LRO can not exist in the GS of such systems, we believe it is true. 
But, if Sa/Sb + (k + l)/(k + 2), the AF and F LROs coexist in the GS. 
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Figure Captions: 

Fig. 1: (a) the ID Hcisenberg AF chain, (b) the QOD Hciscnbcrg AF chain. There arc I sites between site Pi 
and P 2 , and I = 2k + 1, here k = 0, 1, 2.... 

Fig. 2: Spin-spin correlations between site 2 and others (see Fig. 1(a)) for these cases: no side spin (square); 
single side spin at site 2 (cross +), 3 (circle) and 4 (cross x) respectively; two side spins at sites 2 and 4 (triangle); 
two side spins at sites 2 and 3 (diamond). Here, R = Ri — R 2 - The length of chain is 23. 

Fig. 3: Spin-spin correlations between site 2 and others (see Fig. 1(a)) for these cases: bare chain of 21 sites 
(square); 4 side spins at sites 2, 8, 14 and 20 on a chain of 21 sites (diamond); 5 side spins by site 2, 6, 10, 14 and 18 
on a chain of 19 sites (triangle); 8 side spins by site 2, 4, 6, 8, 10, 12, 14 and 16 on a chain of 17 sites (solid circle). 
Here, R = R t - R 2 . 

Fig. 4: p vs. r\. The solid line is obtained by fitting p in linear way. 
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